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ABSTRACT 

Neumayer et al. established the existence of a blue-shifted cloud in the core 
of Centaurus A, within a few parsecs of the nucleus and close to the radio jet. 
We propose that the cloud has been impacted by the jet, and that it is in the 
foreground of the jet, accounting for its blue-shifted emission on the Southern side 
of the nucleus. We consider both shock excitation and photoionization models 
for the excitation of the cloud. Shock models do not account for the [SiVI] and 
[CaVIII] emission line fluxes. However, X-ray observations indicate a source of 
ionizing photons in the core of Centaurus A; photoionization by the inferred flux 
incident on the cloud can account for the fluxes in these lines relative to Brackett- 
7. The power-law slope of the ionizing continuum matches that inferred from 
synchrotron models of the X-rays. The logarithm of the ionization parameter is 
-1.9, typical of that in Seyfert galaxies and consistent with the value proposed 

. The model cloud density depends 
upon the Lorentz factor of the blazar and the inclination of our line of sight 
to the jet axis. For acute inclinations, the inferred density is consistent with 
expected cloud densities. However, for moderate inclinations of the jet to the 
line of sight, high Lorentz factors imply cloud densities in excess of 10^ cm~^ and 
very low filling factors, suggesting that models of the gamma ray emission should 
incorporate jet Lorentz factors < 5. 



for dusty ionized plasmas (Dopita et al.||2002 
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Subject headings: black hole physics - galaxies: active - galaxies: individual 
(Centaurus A) - galaxies: jets - line: formation - relativistic processes 



Introduction 



Neumayer et al. (2007) presented high resolution (0.12") infrared spectral images of 
the central parsecs of Centaurus A, which they utilized to derive important dynamical con- 
straints on gas motions in this region, as well as the mass of the central black hole (see 



Neumayer (2010) for a review of the black hole mass measurements). A feature of the data, 
which has not yet been exploited is the existence of emission line gas clouds approximately 
aligned with the parsec-scale radio jet. Questions which naturally arise in this context are: 
What is the excitation mechanism for this gas and what is the cause of the ~ 100 km s^^ 
velocities of this gas relative to the nucleus? Two candidates for the excitation mechanism 
are: (1) Shock waves driven by the jet into clouds adjacent to the jet and (2) Photoion- 
ization by either emission from an accretion disk or from the jet close to the back hole. 
Shock excitation could arise naturally in the process of entrainment of clouds into the jet 
or simply by jet deflection off these clouds. Considering the possibility of photoionization. 



Lenain et al. (2008) modeled the radio through X-ray and gamma-ray emission from the 
nucleus of Centaurus A in terms of either synchrotron plus synchrotron self-Compton (SSC) 
emission or two component synchrotron emission from a relativistic jet. In a subsequent 
paper they incorporated the newly discovered very high energy (VHE) gamma-ray emission 
from Centaurus A into a revised model in which the 1 — 10^ eV emission is due to synchrotron 
radiation. The significance of a synchrotron explanation for the 1 — 10^ eV emission is that 
it may have a high enough flux (depending upon beaming effects) to be an important source 
of ionizing photons. If an SSC model is used for the 1 — 10^ eV flux, then the flux density 



of ionizing photons is much lower at around 10 eV (cf Lenain et al. (2008)). If one or both 
of the jets in Centaurus A are indeed responsible for the photoionization of the nuclear gas 
clouds then modelling of the emission may lead to interesting constraints on the Lorentz 
factor, which strongly affects the number density of ionizing photons through relativistic 
beaming. 

In this paper, we summarize the observational data in § [2| consider both shock and 
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photoionization models in § [3] and discuss our results in § |4} In an appendix, we derive 
expressions for the ionizing photon density incident on a nuclear cloud in terms of the 
observed flux density and the jet parameters relating to the beaming in both the observer 
and cloud directions. 




Fig. 1. — This figure is reproduced from Figure 7 of Neumayer et al. (2007). The left and 



middle panels compare the [SiVI] and H2 velocity fields. The right hand panel shows the 
differenced velocity field of [SiVI] (i.e. v{[SiVI]) — ^(112) highlighting the blue-shifted cloud, 
coincident with the radio knot SJl ( Hardcastle| 2003). (A color version of this figure is 
available in the online journal.) 



Figure [T| reproduced from Figure 7 of Neumayer et al. (2007) highlights the observational 
data, which are the focus of this paper. The left panel shows the [SiVI] emission spread 
along the inner few arcseconds of both jet and counter jet. The middle panel shows the 



H2 emission in the same region, whose velocity Neumayer et al. (2007) show is indicative of 



disk rotation. Hence the differenced velocity field of [SiVI] in the right hand panel shows 
the non-rotational component of [SiVI]. The striking feature of the [SiVI] emission is the 
feature displaced approximately south-west of the nucleus in right ascension and declination 
by (0.7", —0.6") with a blueshift of approximately 100 km s~^. This feature is referred to 
as the "blue cloud" in the following. A blueshift for this cloud is at first surprising: Most 
analyses of Centaurus A assume that the Southern jet is moving away from us, so that one 
expects a cloud, which may be entrained by the jet, should be redshifted. [Neumayer et al. 



(2007) noted an alternative possibility that the blue cloud could be falling in toward the 



nucleus. In this paper we adopt the view that the blue cloud is on the same side of the jet 
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as us and that it is pushed to one side by the jet - hence the blueshift. The geometrical 
configuration that we envisage is represented in Figure [2] 

Such an interpretation for the cloud velocity raises the question of why the cloud would 
not acquire a substantial forward momentum along the jet direction, leading to a redshift. 



Simulations of jets interacting with inhomogeneous media both in the form of disks Suther- 
land & Bicknell (2007) and spherically distributed clouds Wagner & Bicknell ( 2011a|[b ) show 
that jet-cloud interactions are complex and are often affected by the backflow of radio plasma, 
either from the head of the radio source or from the downstream interaction of the jet with 
other clouds. The back flow impedes forward motion of the cloud so that its resultant motion 
is primarily perpendicular to the jet. 



3. EMISSION LINE MODELS 

In this section we discuss the two principal possible excitation mechanisms of gas in the 
nucleus of Centaurus A: Shock excitation and photoionization. In doing so, we concentrate 
on the blue-shifted cloud (see Figure [T]) since this is separated in velocity from the rest of the 
gas in the nuclear region and offers the prospect of being described by a well-deflned set of 
parameters of number density and photon density. However, we also consider the integrated 
emission from the nucleus. 

The K-band infrared spectra do not offer a large number of emission lines for comparison 
with the models. The principal lines that are evident in the spectra are those of [SiVI] at 
1.9602 /im, Br7 at 2.1656 fim and [CaVIII] at 2.322 fim so that the models are not uniquely 
constrained by the data. However, [CaVIII] is an important discriminant of shock-excitation 
and photoionization and it is possible to derive reasonable conclusions on the relative merits 
of these two mechanisms and also to estimate the cloud densities. 



3.1. Shock Models 

Since we are suggesting that the velocity of the blue cloud is the result of being pushed 
to one side by the passage of the jet it is natural to investigate whether the emission from 
the cloud could be the result of shock excitation. 

The emission from an astrophysical shock wave comprises contributions from the shocked 
gas and the precursor region photoionized by radiation from the shock. Using the MAP- 
PINGS HI emission line code, we have computed a comprehensive grid of shock + precursor 
models for velocities between 160 and 1000 kms~^. The results for [CaVIII] and [SiVI] are 



- 5 - 




Fig. 2. — The proposed geometry of the inner Centaurus A jets, the radio knot SJl and the 
blue cloud. (A color version of this figure is available in the online journal.) 

shown in Figure [3] wherein the shock and precursor regions are shown both separately and 
combined. The effect of varying the pre-shock magnetic field is also shown in Figure [3j 

The compelling message from the shock models is that there is no combination of shock 
velocity and magnetic field, which adequately reproduces both the observed [CaVIII]/Br7 
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Log[ SiVI/Bry ] 

Fig. 3. — Shock models of the ratio of [CaVIII] to Bvy plotted against the ratio of [SiVI] 
to Br7 for the nuclear region of Centaurus A. Combined shock + precursor models as well 
as shock-only models and precursor-only models are shown together. The magnetic field is 
parameterized by /3, the ratio of thermal pressure to magnetic pressure in the pre-shock gas. 
Note that the parameter (3 does not affect the precursor. (A color version of this figure is 
available in the on-line version.) 

and [SiVI]/Br7 line ratios. The adjacency of gas to the jet combined with its velocity offsets 
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with respect to the core indicates that its kinematics may be affected by the jet. However, 
it does not appear that shock; waves associated with a jet-cloud interaction dominate the 
emission. 



3.2. Photoionization Models 



We now consider photoionization models for the emission from the nuclear region of 



Centaurus A. The point of view that we take here, is inspired by that proposed by [Pedlar 
et al. (1989) and Taylor et al. (1992) for Seyfert galaxies. Clouds of moderate density 



are impacted by the radio jets producing higher density radiatively shocked regions. The 
shocked gas is photoionized by the nuclear continuum, producing the observed emission line 
spectrum. 



3.2.1. Photoionizing Flux 



One possible source of the ionizing continuum is the nuclear sub-parsec-scale jet. In 
order to assess the implications of such a source, which may be relativistically beamed, we 
need to determine the flux of ionizing radiation incident on the cloud given the radiation 
that we observe along a line of sight, which is inclined at a different angle to the velocity of 
the radiating plasma. 

The high energy emission of the Northern jet from hard X-rays through to very high 



energy (Tev) gamma rays has been modeled by Lenain et al. (2009 ) in terms of a Synchrotron- 



Self-Compton (SSC) model. This was an update of an earlier model (Lenain et al. 2008) in 



which both synchrotron and inverse Compton models were considered for the X-ray emission 



alone. However, the detection of VHE 7-ray emission from Centaurus A (Aharonian et al. 



2009) has made the synchrotron interpretation more appealing. This is important since the 



inverse Compton model has a much lower ionizing flux. In the model that we develop here, 
we assume that the base of the Southern jet has a similar spectrum to the Northern jet and 
that this provides the flux of photoionizing radiation. 



Figure 1 in Lenain et al. (2009) shows that the continuum region that is of interest 



here, namely the ionizing continuum with energies in excess of approximately 10 eV in the 
rest frame of the irradiated gas, is only constrained by the observational data in the range 
~ 1 — 100 keV. Hence we rely on the modeled spectrum in the region ~ 1 — 10^ eV. The 
model spectral energy distribution (SED) between approximately 1 eV and 100 keV (the 
power-law synchrotron component indicated on their Figure 1) has a flux spectral index 



- 8 - 



of a ~ 0.39 {Fy oc v °). This is slightly flatter than usual (say ~ 0.6 — 0.7) for a non- 
thermal spectrum. Nevertheless, values of 0.4 and 0.3 have been inferred for the synchrotron 



component of the blazar emission in MKN 501 by Bicknell et al. (2001) and Konopelko et al. 



(2003) respectively. 



Lenain et al. (2009) present three different models. In the following we use the model 



which most accurately reproduces the X-ray data, the thin line in their Figure 1. 



Lenain et al. (2009) infer a Lorentz factor of 15 for the Centaurus A jet. Lorentz factors 



of this magnitude are controversial in view of the fact that Lorentz factors estimated on 
the parsec scale from radio observations are always lower than this (e.g. Finer et al. 2008) 



and Centaurus A is no exception: The pattern speeds of knots in the Centaurus A jets are 
of order 0.1 c (Tingay et al. 1998) - much less than the superluminal pattern speeds seen 



in many quasars, although Tingay et al. note that there is evidence for a faster underlying 
flow. Estimated high Lorentz factors in blazars may point to significant deceleration between 
the sub-parsec and parsec scales or to localized fast-moving regions of the flow and these 
issues are currently unresolved. However, the inference of high Lorentz factors and the 
observation of superluminal motions in AGN in general, means that we need to take account 
of relativistic beaming. As noted above, in Centaurus A the radiation field observed by us 
should be different from the radiation field intercepted by the blue-shifted cloud and it is 
important to allow for this in the photoionization models. 

The various relationships required to relate the observed flux density to the ionizing 
photon density at the cloud are derived in the appendix. We summarize the main results 
here: 

Our model is summarized in Figure |2} We envisage the inner jet, within say about 
100 gravitational radii (~ 20 mpc) as the source of ionizing photons (the inner jet); let 
Del be the distance of the cloud from the source; let be the angular diameter distance 
of Centaurus A; let the Doppler factors of the radiation received by observer and cloud 
respectively be (5obs and 5ci and let z ^ 0.0018 ( Graham|1978 ) be the redshift of Centaurus A. 
If v is the frequency of a photon intercepted by the cloud, it is emitted by the jet with a rest 
frame frequency S'-^u and a photon emitted with this frequency reaches the observer with a 
frequency Vohs = 5obs(l + ^)"^ 
jet as a function of the observing frequency i^obs- 



6^^ u. Let -Fobs(t'obs) be the observed flux density of the inner 



With these definitions, the photon density per unit frequency at the cloud (derived in 
the Appendix) is 



Da 
Del 



Sclil + z) 



obs 



-Fobs('^obs(l + Z) 



(1) 
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Let z/q = 3.28 x 10^^ Hz be the frequency corresponding to the Rydberg hmit. The total 
number density of ionizing photons at the cloud (also derived in the appendix) is 



n 



ph 



nph(z/) diy 
1 [Da 



5M + z) 



1 3 



5, 



obs 



<Sobs(l + ^)~^5c/^0 



-^obs(l^obs) 



dv, 



obs 



(2) 



Between about 10^^'^ and an upper frequency z/„ ~ 10^^'^ Hz the modeled spectrum may 
be described as a power-law F(z/obs) ~ -^(t'o)(^'obs/ t'o) " with F(z/o) 
and a ~ 0.39. Hence, 



1.35xlO-27ergss-iHz-i 



obs 



>o) 





2 

) 




3+a 




'5cl(l + Z) 








l^obs 




1- 


^obs 


1^0 



(3) 



Since z/^j 3> z/q, the contribution of the second term in brackets is minor. 



Let the inclination of our line of sight to the jet be ^obs, the angle between the direction 
of a ray from the core through the center of the cloud and the jet be ifj^x and the (projected) 
angular displacement of the cloud from the core be ^ci, the geometry of the model (see 
Figure |2|) implies that the ratio of angular diameter distance to cloud distance is given by 



Da _ sin(6'obs + ^ci) 



cl 



(4) 



(This equation is also derived in the Appendix.) Furthermore, the ratio of Doppler factors 
is given by: 

5c\ _ 1 - /3 COS gpbs 
5ohs 1 - /3 COS Ipcl 

As expected, the expression for the ionizing photon density has a strong dependence on 
the ratio of Doppler factors, when ifjcx is small. We estimate ipci assuming that the cloud 
is adjacent to the jet and that is depth along the line of sight is the same as its transverse 
diameter 



0.35". The expression for ifjci is derived in the Appendix (see equation A18). 



3.2.2. Isochoric Photoionization Models 

In framing models for cloud photoionization one question to resolve is whether an iso- 
baric or isochoric model is more appropriate. Here, an isobaric model is relevant if the 
expansion time scale for the gas heated by photoionization is less than the time scale for 



- 10 - 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 




Log U = +0.0 : 



- a = 0.0 - 
I I I I I I I I I I I I I I I I I I I I I I I I 

-0.5 0.0 0.5 1.0 

Log[SiVI/Bry] 

Fig. 4.— The line ratios [CaVIII]/Br7 and [SiVI]/Br7 resulting from a grid of MAPPINGS 
isochoric, photoionization models for solar metallicity. The spectral indices are marked on 
the different curves and the dots on the curves represent increments of 0.1 in the logarithm 
of the ionization parameter. Data points corresponding to (a) the blue cloud and (b) the 
total nuclear emission are also marked. 
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cloud shredding resulting from the jet-cloud interaction. The photoionized cloud would ex- 
pand at approximately the sound speed corresponding to a temperature of 10^ K, that is, 
about 12kms~^. If we take the anomalous velocity of the blue cloud ~ lOOkms"^ as indicative 
of shocks driven into it by the jet, then the shock-shredding time scale ~ 18pc/100kms~^ ~ 
1.8 X lO^yrs is approximately a factor of eight shorter than the expansion time scale. Hence, 
we mainly consider isochoric, constant density, models in the following. 

We have calculated a grid of models for solar metallicity using the solar abundance scale 



of Grevesse et al. (2010). The models are parameterized by the ionization parameter, U, 
the ratio of the ionizing photon density, riph to the particle density n {U = ripi^/n) and the 
power-law index, a of the photoionizing flux density [F^, oc The models are ionization 

bounded and are for a density n = 10^ cm~^. For this fiducial density, the extent of the 
ionization bounded region ^ 10^^'^cm and fits comfortably within the projected major axis 
]^g20.7 ^]^g blue cloud. Below, we consider lower filling factors, involving higher 

densities. However, the model line ratios remain independent of density for all the densities 
considered here. 

Figure |4] shows the results of the grid of models for the line ratios [SiVI]/Br7 and 
[CaVIH]/Br7. The best-fit model is represented by (logf/, a) ~ (—1.9,0.42). The best- 
fit value of a is in practice indistinguishable from the value inferred from the synchrotron 
emission fit to the X-ray data, strongly supporting the notion of photoionization by jet 
emission. 

We note that the best fit ionization parameter of logf/ ~ —2.0 is in rough agreement 



with the self-limiting radiation pressure dominated dusty AGN nebula models of Dopita et al. 



( |2002D ; [Groves et al] ( |2004a||bp . In future work it may be fruitful to explore the possibility 



of such dusty models in the present context. However, with a scarcity of diagnostic lines such 
as [OIII] A5007, and a lack of constraints on the properties of the dust, and any depletion 
factors that may be present in the blue cloud (affecting the abundances of both Silicon and 
Calcium), additional model parameters are not justified by the present spectral data. 

If more data on the cloud composition, dust properties and more spectral line mea- 
surements were available, then radiation pressure dominated dusty models may improve the 
models and place better constraints on the Br7 emission for example. This would affect the 
filling factor of the cloud (see below). 

In our isochoric models, the ratios of [SiVI] and [CaVIII] to Br7 are primarily determined 
by the slope of the ionizing continuum and the level of the emission is determined by the 
density. Hence, we can estimate the cloud number density, n from the ionization parameter, 
U, and the ionizing photon density, riph (using n = U^^riph). Since the ionizing photon 
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density depends upon the observed continuum, the angle of inchnation of the radio jet and 
its Lorentz factor, through equation (|3]), we can assess the effect of these quantities on the 
inferred density of the cloud. Note that our estimates of the number density would not be 
substantially revised by dusty models since such models predict a very similar ionization 
parameter to what we have inferred here. 



3.2.3. Filling Factor. 

The filling factor of the Hydrogen-line emitting region of the cloud can be estimated from 
the Br7 luminosity, L(Br7). Let a(Br7) be the effective Case B recombination coefficient 
for Br7, / the volume filling factor, Ue and nu the electron and Hydrogen number densities 
and Vci the cloud volume, then 

L(Br7) = a(Br7) / neUH Ki- (6) 

We estimate a cloud volume 

Kl -D\^ (7r/6) ^n^ajCin/ Sm ^obs (7) 

where ^maj ~ 1-1" and 6*111111 ~ 0.7" are the angular sizes of the major and minor axes of 
the cloud respectively; we assume that the depth of the cloud is similar to the length of the 
minor axis. 

In using Equations (|6| and ([T]) to estimate the filling factor, we use the measured 
value of the Br7 fiux of 2.42 x 10~^^ergs s~^ cm~^ and adopt a value of 2.67 x 10~^^ for 
the Br7 emission coefficient, appropriate for the average temperature of 12, 500 K obtained 
in the photoionization models. Note that the estimate for the filling factor depends upon 
the assumption that the cloud depth is similar to its width. If the depth were an order of 
magnitude smaller, for example as a result of ablation by the jet, then the filling factor would 
increase by an order of magnitude. 

Our estimates of the density of the photoionized gas and the related filling factor both 
as a function of Lorentz factor and ^obs are presented in the left and right panels of Figure [5] 
For all 6'obs, the density increases quite rapidly with Lorentz factor. For 6^obs = 40°, 60° and 
80°, the density increases to values in excess of 10'^ cm~^. In keeping with this, the filling 
factor decreases quite rapidly with Lorentz factor for all inclination angles but especially for 
6'obs > 40°. The strong dependence of the density and filling factor on Lorentz factor results 
from the effects of beaming and the high inferred photon density at the cloud when the line 
of sight to the observer is greater than approximately 30° to the jet axis (see equation ([3])). 
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Fig. 5. — Left panel: Estimate of the cloud density from the ionization parameter and the 
ionizing photon density, which is a function of the jet Lorentz factor and the inclination of the 
jet to the line of sight. Right panel: Estimate of the filling factor from the Br7 luminosity. 



That is, for moderate values of the angle of the observer to the line of sight, the observer 
sees a considerably lower intensity than that which is incident on the cloud. 

Various inclination angles have been proposed for the Centaurus A jet but most are 
within the range 50° < 6'obs < 80° proposed by Tingay et al. (1998). The exception is 



Hardcastle et al. (2003) who estimate an angle of inclination of approximately 15° if the 
inner Very Large Array (VLA) jets are intrinsically symmetric. However, they note the 
inconsistency of this estimate with other estimates of inclination, leaving open the question 
of the asymmetry of the inner VLA structure. 

What is the evidence for gas with densities in excess of 10^ cm~^? Using a spherical 



free-free absorption model for the very long baseline interferometry (VLBI) jets, Tingay 
& Murphy (2001) derived a number density ~ 9 x 10^ T^ '^^ where the temperature of the 



absorbing gas is 10 T4 K and the radius of the free-free absorbing sphere is 0.016 pc. Wild 



et al. (1997) determined gas densities ^ 1 — 3 x 10^ cm ^ from CO observations of the 



Centaurus A dust lane, at a spatial resolution of 44". Subsequent work (Wild & Eckart 



2000) using HCN and CS observations at a resolution of 54" showed that approximately 
10% of the molecular gas in the nuclear region has densities in excess of 10^ cm~^. Thus, 
there is good evidence for dense gas clouds in the nuclear region of Centaurus A, associated 
with the dust lane. Future ALMA observations should provide more information on the 
spatial distribution of these clouds. However, we do not expect the cloud densities at a few 
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parsecs to be higher than the 10 cm inferred by Tingay & Murphy (2001) near the core 



so that we adopt a conservative upper estimate on the density of 10 cm . 

It is only for inchnations < 50° that the constraint n < 10^ cm~^ is reahzed for all 
Lorentz factors up to the estimate of F = 15 by Lenain et al. (2009). This density constraint 
is consistent with the Lenain et al. (2009) models in which a low inclination angle, 6'obs = 25°, 



is assumed but not with the range of inclinations discussed above. 

Let us now consider the dependence of filling factor on Lorentz factor. While describing 
the density structure of a photoionized region using this parameter is not ideal, we adopt if 
here for the purpose of approximate comparisons. In the narrow line region of active galaxies 
the value inferred for the filling factor is typically ~ 10^^ — 10~^. Pedlar et al. (1989) and 



Taylor et al. ( 1992 ) explained this low value in terms of the dynamics of radio plasma driven 



radiative shocks, which produce high density, but low volume filling factor, shocked regions. 
The high density regions are photoionized by the nuclear radiation and the Pedlar et al. and 
Dyson et al. models explain why photoionization models require a low filling factor. Pedlar 
et al. and Dyson et al. formulated their models to explain the observed properties of Seyfert 
galaxies, and the physical situation envisaged here is similar to their Seyfert model: Gas 
shocked by the jet is photoionized by the central source producing regions of photoionized 
gas with a low filling factor. Hence a comparison of the filling factors in the narrow line 
regions of Seyfert galaxies and those deduced here for Centaurus A is useful. We therefore 
adopt a conservative lower limit on the filling factor of approximately 10~^. 

When we examine the dependence of filling factor on Lorentz factor (right panel of 
Figure |5| we see that the inferred filling factor decreases rapidly with Lorentz factor. This 
results from the increase of the beamed photon density and atomic density with Lorentz 
factor. The effect is highlighted for moderate inclinations. For acute inclinations the value 
of the filling factor is acceptable up to F 10. However, for moderate inclinations, filling 
factors < 10~^ are only possible for Lorentz factors less than a few. 



4. Conclusions and discussion 

The discovery of blueshifted line-emitting gas on the side of the nucleus that we would 
normally associate with redshifted outflow has led us to consider detailed models for the 
location, kinematics and excitation of this gas. 

The features of the model that we have investigated in detail include the following: 
(i) The blue-shifted cloud has been impacted by the Southern jet and is located close to 
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the jet and between the jet and us. 



(ii) The gas has been shocked by the jet producing dense sub-clouds/filaments of gas but 
the primary excitation mechanism for the hne emission is not shock excitation but 
photoionization by uhraviolet and X-ray emission from the core of Centaurus A. 

(iii) The X-ray emission from the core combined with the observations and models of the 



VI 



vn 



high energy gamma-rays ( Lenain et ah 2009 ) suggest that the ionizing flux originates 



from the Southern jet. We have assumed that the Southern jet flux is similar to that 
of the Northern jet. 

(iv) Photoionization models of the emission from the blue cloud require a power-law slope 



almost identical to the Lenain et al. (2009) models of the UV-X-ray emission from the 



core of Centaurus A, in order to reproduce the [CaVIII] and [SiVI] fluxes relative to 
Br7. 

The power-law slope is characteristic of the high-energy slopes of core jets in blazars. 

The last two results indicate very strongly that ionizing emission from a relativistic jet 
in Centaurus A is responsible for the excitation of not only the blue cloud but other 
nuclear line emission in Centaurus A. 

Consideration of both the number density and volume filling factor implied by the 
photoionization models and the flux of Br7, strongly suggest, that the Lorentz factor 
of the jets in Centaurus A is lower than indicated by the high energy models. This is 
not surprising since it has been an unresolved problem for some time that blazar models 
usually imply higher Lorentz factors than those inferred from VLBI observations on 
the parsec scale. 



As we have seen in § 3.2 precise conclusions concerning the Lorentz factor of the jet in 



Centaurus A (and consequently the viability of models for the high energy emission) depend 
strongly on our viewing angle and the interpretation of the milliarcsecond and arcsecond 
radio data. Broadly speaking, if the inclination of our line of sight to the jet axis is as acute 
as the interpretation of the arcsecond radio data may suggest (Hardcastle et al. 2003[ ) then 
our photoionization models are consistent with a reasonably high Lorentz factor (for both 
jets). If, on the other hand, the inclination is in the range of 50° — 80° inferred from the 
millarcsecond radio data (Tingay et al. 1998), then low to moderate Lorentz factors (F < 5) 



are implied . The is consistent with the lack of evidence for very high Lorentz factors in the 



form of superluminal proper motions (Tingay et al. 1998) such as are found in quasar jets 



for example (Kellermann et al. 2004). 
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The physical distance, of the blue-shifted cloud from the nucleus is determined by 
the model parameters, specifically 6'obs; the angle between the observer and the Northern 
jet direction, and ipci, the angle between the direction of the southern jet and the direction 
of the cloud from the nucleus, together with the projected distance from the nucleus ~ 
16.5 pc corresponding to an angular separation of approximately 1 arcsec. (See Figure |2] 
and Equation A15 ) For the most acute angle of inclination we have considered, ^obs = 20°, 
-Del ~ 36 pc. For the range of angles referred to in the above paragraph 50° < 6'obs < 80° the 
distance of the cloud from the core varies very little: 18 pc > Dd > 17 pc. Note however, 
it is only in the case of the photoionization models that these parameters are relevant. The 
main parameter in the shock models is the shock velocity. 

Our conclusion that the blue cloud is photoionized by beamed emission from the core 
invites the question as to the relevance of photoionization for other clouds on the kiloparsec 



scales in the interstellar medium of Centaurus A (Blanco et al. 1975 Peterson et al. 1975 



Graham & Price 1981; Morganti et al. 1991). In particular Morganti et al. (1991) argued 



that the region referred to as the inner filaments are photoionized by beamed radiation 



from the core. On the other hand, Sutherland et al. (1993), presented a model in which 
the filaments are locally shock-excited through interaction with the radio plasma. More 



recently Sharp & Bland-Hawthorn (2010) have found an ionization cone in the inner galaxy 
of Centaurus A aligned with the radio jet as well as other nearby emission line regions, which 
are excited by star formation. In the light of the calculations presented in this paper, we can 
consider some examples based on typical parameters which indicate the conditions under 
which photoionization may be feasible for the different regions well outside the core. Note 
that these regions are all to the North of the nucleus so that we do not need to assume the 
equivalence of Southern and Northern jet ionizing fluxes. 



Region 


eel 


r 


^cl 


n 


u 






arcsec 




degrees 


cm~^ 






Ionization 


25 


2 


17.5 


100 


4.1 X 10" 




Cone 


25 


4 


17.5 


100 


4.9 X 10- 


-3 




25 


5 


17.5 


100 


7.9 X 10- 


-3 


Inner 


550 


2 


3.5 


30 


2.2 X 10- 


-5 


Filaments 


550 


4 


3.5 


30 


1.8 X 10- 


-3 




550 


5 


3.5 


30 


6.8 X 10" 


-3 


Outer 


910 


2 


10.4 


10 


5.9 X 10- 


-6 


Filaments 


910 


4 


10.4 


10 


2.0 X 10" 


-4 




910 


5 


10.4 


10 


4.8 X 10" 


-4 



Table 1: Estimates of Ionization Parameter in Different Regions of Centaurus A 
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In Table [T] we present some indicative calculations of ionization parameters for these 
different regions for various values of jet Lorentz factor, nominal values of number density 
and a jet inclination of 60°. The angle between jet and cloud directions were estimated as 



follows. For the ionization cone Sharp & Bland-Hawthorn (2010) estimated a half-opening 



angle of 20°; deprojected this is 17.5°; our estimate of ionization parameter refers to the edge 
of the cone. The inner and outer filaments have position angle differences with respect to 
the milliarcsecond jet of 4° and 12° respectively. Since we are only determining indicative 
numbers, we assume that the inclination of the cloud direction to the line of sight is the 
same as the inclination of the jet to the line of sight (i.e. ^ci = ^obs) and we calculate the 



angle between the jet and the cloud direction, ipci, using equation (A18) 



The cloud densities assumed for the inner filaments are the same as the estimates of 



Morganti et al. ( 1991 ); the other densities are nominal but the resulting ionization parameters 



can easily be scaled for other densities. 

For the spectra that are observed in Centaurus A the required ionization parameters 
are of order a few xlO^^ to 10^^. Indeed Morganti et al. (1991) estimated a reference value 
from their variable density models for the inner filaments of 5 x 10^^. Our estimates for 
the ionization cone and the inner filaments are in this range for jet Lorentz factors of four 
and five but not for the lower Lorentz factor of two. If the outer filaments are to fall in this 
range the density would need to be an order of magnitude lower. Hence the estimates of jet 
parameters derived from models of the blue cloud are consistent with photoionization of the 



Sharp & Bland-Hawthorn (2010) ionization cone. However, the situation is not as clear for 



the outer filaments. 
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supported by the Australian Research Council Discovery Project DP0664434. NN acknowl- 
edges the support by the DFG cluster of excellence Origin and Structure of the Universe. 
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A. RELATIONSHIP BETWEEN OBSERVED EMISSION AND BEAMED 

PHOTOIONIZING FLUX 

We often see high energy X-ray and 7-ray emission from the cores of active galaxies, 
which is interpreted as emission from a relativistic jet. Potentially, this emission can ionize 
clouds along the direction of the jet where it is more highly beamed. However, as a result 
of the beaming pattern of the emission, such clouds see a different intensity to that at the 
observer. In this Appendix we provide a ready way of estimating the ionizing photon density 
given parameters such as the Lorentz factor of the jet, the directions of the observer and 
cloud and the respective distances of cloud and observer. 



This preprint was prepared with the A AS IM^jX macros v5.2. 
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A.l. Ionizing Photon Density 

We use the following symbols and definitions: 

(i) I{i>) is the specific intensity of a ray propagating from a region (blob) in the jet through 
the cloud. 

z/q = 3.28 X 10^^ Hz is the frequency of the Lyman limit 

dfl is the elementary solid angle of rays emanating from the relativistic blob and 
intersecting the cloud 

The rest frame of the blob is referred to with a prime. 
j'(z/') is the (isotropic) emissivity in the blob rest frame. 
x' is the path length along a ray in the blob rest frame. 

dA is the element of area, in the cloud frame, of the surface S bounding the relativis- 
tically moving plasma; dA' and E' are the corresponding quantities in the blob rest 



frame. As spelled out in Lind & Blandford (1985) S and S' are congruent; S is E' 



rotated by relativistic aberration and dA = dA'. 

dQ' is the elementary solid angle in the blob rest frame. 

(3 = v/c for the moving plasma; F = (1 — /3^)~^/^ is the Lorentz factor. 

ipc\ is the angle of the ray passing through the center of the cloud measured with respect 
to the direction of the jet axis. 

6ci = r^^(l — f3cosipci)~^ is the Doppler factor of the emitting jet plasma as viewed 
from the cloud. 

The frequency, z/, in the cloud frame and the frequency, u' in the blob rest frame, are 
related hy u = 6ci v' ■ 

Del is the distance of the cloud from the relativistic blob. We assume that D^i is large 
compared to the dimensions of the emitting region. 

dV is the elementary volume of the relativistic blob. 

Da is the angular diameter distance of the galaxy containing the relativistic jet. 
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(xvi) 5ohs — r ^(1 — ;5cos^obs) ^ is the Doppler factor of the source along a ray to the 
observer; ^obs is the direction of the observer with respect to the blob's velocity. 

(xvii) The frequency u' of a photon emitted in the jet rest frame and the observed frequency 
i^obs are related by i^obs — <^obs(l + where z is the redshift of the galaxy. 



The number density of ionizing photons (per unit frequency and total) are, respectively: 



hp 



hv 



dn 



dv 



(Al) 
(A2) 



We relate these quantities to the parameters of the relativistically moving emission region 
as follows. 



The intensity in the cloud frame and jet rest frame are related by: 
where 



(A3) 



(A4) 



Ray 



and the integral is along a ray through the relativistically moving emitting region TZ. That 
is, 

j\5;,^v)dx' (A5) 



m = Si 



Ray 



Hence, the number density of photons per unit frequency is 



7 ^w,,; 

./Ray hu 



dfl 



(A6) 



The number density of ionizing photons in the cloud frame may be determined by 
integration over the volume of the relativistic emitting region TZ. Since dfl = dA/D^ = 
dA'/Dl, then dx'dVL = dx'dA'/D^ = dV/D^ and 



ripM = - 7^ 



hv 



dV. 



(A7) 



The integral over solid angle and path length has been replaced by an integral over comoving 
volume. We integrate over frequency to obtain the total number of ionizing photons. Thus, 



n 



ph 



Blob 



hv 



dV 



dv 



(A8) 
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A.2. Flux Density at Observer 



For completeness, and in order to relate observed flux density to ionizing photon density 



we repeat the analysis in Lind & Blandford (1985) in the following expression for the flux 



density. This enables us to estimate the volume integral of the emissivity which appears in 



the expressions for the ionizing photon density, equations (A7) and (A8). 



The flux density of the blob, at frequency i^obs, as measured by the observer is: 



Fiv, 



obs / 



-^obs(i^obs) d£L 



(A9) 



As in previous analyses (e.g. Lind & Blandford 1985; Begelman et al. 1984) we determine 
the observed flux density in two stages, from blob to a point well outside the blob in the 
galaxy rest frame and from galaxy rest frame to observer. The flrst stage involves relativistic 
effects; the second stage cosmological effects. The result for the observed flux density is: 



F(z/, 



obs ) 



1 



6, 



obs 



where is the emissivity in the jet rest frame. 



Oobs 



(AlO) 



A. 3. Using the Observed Flux Density to Estimate the Ionizing Photon 

Density 

We now use the observed flux to determine the volume integral of the emissivity and then 
the ionizing photon density. The combination 5obs(l+2;)~^5^^z^ which appears in the following 
reflects the fact that a photon which intersects the cloud with frequency z/, originates from 
the blob rest frame with frequency Photons emitted with this frequency reach the 

observer with frequency 5obs(l + z)~^5~^v. 



From equation (AlO) for the flux density, we have: 



When we insert this expression into equation (A7) for the number density of ionizing 
photons per unit frequency we obtain: 

c \Ddl V c)obs / hu 
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and the number density of ionizing photons is: 

We now change the integration frequency to the observer's frequency z^obs = 5obs(l + ^)^^5ci^^ 
to obtain: 

1 /daV f{i+z)5,iY r ^Kbs) 



'^ph = - 7f , / -ir^ ^'^obs (AM) 

C \Uc\J \ Cobs / ^5obs(l+^)-i5e^^t'o '^'^obs 

We now estimate the factor D^/ Dd appearing in this equation. Let the angular separation of 
the cloud from the nucleus in the plane of the sky be ^ci and the projected linear separation 
be -Dci,p- Then, referring to Figure |2] for the source geometry, the polar angle of the ray 
intersecting the center of the cloud is 6'ci = tt — 6'obs — V'ci and the projected cloud distance 
from the core is 

^ci,p = Da^cX = Del sin6'ci = D^i sin(6'obs + V^ci) (A15) 

Hence, 

Da sin(6'obs + V'ci) 



Del ^cl 

Thus, the number density of ionizing photons at the cloud is: 



(A16) 



lsm\e^i,, + ^ei) f{l + z)6eiY r ^Kb: 



riph = T2 — 7 / —r c^'^obs (A17) 

This form makes it straightforward to use the observed ionizing spectrum to evaluate 
the photon density at the cloud as a function of the jet Lorentz factor and the angular 
direction of the rays intersecting the cloud. 

The angular parameter ipd depends upon the angle of inchnation 6'obs- Consider a 
cloud which is adjacent to the jet and assume that it is approximately spherical. The angle 
between a ray through the cloud center and the jet axis is determined by the sum of the 
jet and cloud radii and the distance of the cloud from the core. The jet and cloud radii 
are independent of the angle of inclination but the distance from the core is affected by 
projection. Let $ci ~ 0.35" and $jct ~ 0.35" be the observed angular diameter of the cloud 
and jet respectively. Then the angle, ipci between the ray through the center of the cloud 
and the jet axis is given by 



■001 = tan 



sin 9, 



obs 



— sm^obs (A18) 



2^01 

where, as above, ,^ci is the projected angular separation of the cloud from the core. In our 
calculations we round up the value of ipci to the nearest degree. 



